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TECHNICAL NOTE 4088 

PRACTICAL SOLUTION OF PLASTIC DEFORMATION PROBLEMS 
IN ELASTIC -PLASTIC RANGE 
By A. Mendels on and S. S. Manson 


SUMMARY 

A practical method for solving plastic deformation problems in the 
elastic-plastic range is presented. The method is one of successive 
approximations and is illustrated by four examples which include a flat 
plate with temperature distribution across the width, a thin shell with 
axial temperature distribution, a solid cylinder with radial temperature 
distribution, and a rotating disk with radial tenp>erature distribution. 


INTRODUCTION 

The calculation of stresses in structural components in which plas- 
tic flow is considered is currently of great interest in order to take 
full advantage of the load-carrying capacity of available materials . 
Little attention, however, has been directed at providing simple, general 
methods which can be applied by the engineer toward the solution of prac- 
tical problems. This report presents such a method and its applications 
to several problems of current interest. Although use is made of a 
technique arising in the theory of integral equations, no knowledge of 
integral equations is required, and the mathematics involved is well 
within the scope of the practicing engineer. 

The method, which is one of successive integrations, is illustrated 
for four different thermal stress problems which include the flat plate 
with temperature variation along the width, the thin cylindrical shell 
with axial temperature distribution, the solid cylinder with radial tem- 
perature distribution, and the rotating disk with radial temperature dis- 
tribution. The techniques illustrated are not, however, limited to ther- 
mal stress problems. The first three problems considered involve small 
plastic strains, that is, on the same order of magnitude as the elastic 
strains . The fourth problem involves strains on the order of 1 percent . 
The deformation theory of plasticity with the Von Mises yield condition 
is used. Other yield conditions, however, could be used. 
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METHOD 

The determination of stresses and strains in a body consists of com- 
bining the equations of equilibrium and compatibility with the "stress- 
strain" relation and integrating the resulting equations. For both elas- 
tic and plastic problems the same equations result for equilibrium and 
campatibilityj the difference in the two cases consists of the stress- 
strain relation. In the elastic case a linear relation applies between 
stress and strain, thus resulting in linear differential equations which 
can be solved by standard methods . The occurrence of plastic flow great- 
ly complicates the problem by introducing a nonlinear stress-strain rela- 
tion resulting in a nonlinear differential equation for which direct 
solutions may be difficult or impossible to obtain. The present report 
demonstrates a relatively simple method for obtaining approximate solu- 
tions to the resulting nonlinear differential equations whereby the equa- 
tions are first converted to nonlinear integral equations and then solved 
by t he established mathematical technique of successive integrations. 

The principle underlying the successive integration method of the 
solution of plastic flow problems can best be illustrated by an example 
which is treated in greater detail in the section EXAMPLES. For a flat 
plate subject to a temperature variation along the chord 

T = 600(y 2 - l/3) + Tq the total strain e due to stress is given by 

e =• -0.00570(y 2 - l/3) + / e ep ^ (l) 



where the plastic strain e is a function of the total strain e 
(fig. l). (All symbols are defined in appendix A.) 

A convenient method for solving equation (l) is to approach the 
solution in successive steps by organizing the computations so that the 
bothersome nonlinear terms are treated not as unknowns but as known 
quantities determinable from a previous iteration. For example, as a 
zeroth approximation, let it be assumed that e ep = 0 for all values of 

y. Then equation (l) permits direct computation of. e for all values 
of y. Once the total strain e is determined, the plastic component 
of the strain e„_ may be obtained by inspection of the stress-strain 

curve (fig. l) or by simple confutation. It is thus possible in the next 
computation to include the e ep terms as determined from the first com- 
putation, and thus somewhat more accurate values of e can be obtained. 
These in turn lead to more accurate values of e ep , and the process is 

repeated, each time with the plastic strain terms determined in the pre- 
vious iteration being treated as known quantities. When successive iter- 
ations lead to no change in e or e ep , it manifestly makes no difference 
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whether the e terms apply to the current or previous iteration. 

Equation (l) is now satisfied to the desired degree of accuracy, and the 
correct solution is therefore obtained. 


The simple example of the flat plate just illustrated involves uni- 
axial stress, hence it is possible to determine directly plastic strain 
once the total strains are known. In most cases of practical interest 
the stresses are biaxial or triaxial, and the formal procedure for carry- 
ing through the computations can be illustrated by referring to the case 
of the long solid cylinder with radial temperature distribution, which is 
also treated in greater detail in the section Long Solid Cylinder. By 
manipulating the equilibrium and compatibility equations and the stress- 
strain relations, expressions for the total strains in three principal 
directions are derived in integral equation forms of the type 



g(r , e^. , e e , e z , £rp , e 0p )dr 


( 2 ) 


with similar expressions for e 0 and e z . If an attempt is made to write 

a direct relation between the elastic and plastic strains derived from 
plasticity laws and this relation is substituted in equation (2), compli- 
cated nonlinear integral equations result which do not readily admit solu- 
tion. These complications can again be avoided by the process of succes- 
sive integrations. The zeroth approximation for the total strains e r , 

£ q , and e z is determined by neglecting the integral term, thereby e r 
is made equal to f(r) with similar assumptions for e 0 and e z . For 
these values of total strains the plastic strains are determined as will 
be discussed presently. These values of plastic strains are now treated 
as known values and substituted into equation (2) to determine first 
approximations to total, strains e r , e^, and e z from which new approx- 
imations to the plastic strains can be determined. The process is re- 
peated as many times as necessary until successive approximations show 
sufficiently little change in total or plastic strains to permit being 
considered as converged to the correct solution. 

It may be noted that the method of successive integrations is not 
completely new in application to plastic flow problems. Ilyushin’s 
treatment of the thin shell (ref. l) is essenti ally a successive inte- 
gration process similar to that discussed in this report. Although he 
regarded the successive iterations as a series of artificial elastic 
problems, the mechanics of the computations are identical to those re- 
sulting when the problem is treated strictly in its mathematical sense 
of successive integrations of a nonlinear integral equation. In the 
present paper the thin shell problem is treated without the restrictions 
imposed by Ilyushin of linear strain-hardening and of complete incompres- 
sibility in the elastic-plastic range. 
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In order to determine new values for the plastic strains and 

e a from the values of the total strains e . e Q , and e_ as calculated 
& p r 8 z 

by equations such as equation (2), a stress-strain relation for biaxial 
or triaxial stresses in the elastic -plastic range is needed. It will be 
seen that the validity of the method does not depend on the precise form 
of the stress-strain relations, but for illustrative purposes it was nec- 
essary to select specific relations . The relations of the deformation 
theory of plasticity were therefore used. Appendix B shows that by in- 
troducing the concept of equivalent total strain, the plastic strains can 
be computed from the total strains in a simple fashion. Thus, the equiv- 
alent total strain e t may be defined as follows: 


~et 


= “ ^ \J( e x ~ e 0^ 2 + ( e r " e z) 2 + ( e 0 ” e z^ 


Then the plastic strains are given by: 


( 3 ) 


= k IT 2 ( 2( V - e n " O 


"rp 3 e 


et 


et 


( 4 ) 


where is the equivalent plastic strain as shown on the uniaxial 

tensile curve of figure 1. Furthermore, it is shown in appendix B that 
the equivalent total strain e gt can be written as follows: 


_ 2(1 + v) fe 
e et 3 E 


+ e 


ep 


( 5 ) 


where the equivalent stress o 0 is the ordinate of the uniaxial stress- 
strain curve as shown in figure 1. Since, for a given value of e ep> 
cr 0 can be directly determined from this stress-strain curve, equation 
(5) pe rmi ts the direct construction of a curve of € et against e e p as 

shown in figure 2. The computation of the plastic strains thus becomes 
very simple. For a given set of values of total strain, is computed 

from equation (3) . For this value of e gt> e ep figure 2, 

and e-pp and e 0p are then computed from equations (4). These values 
of and €gp are substituted into equations such as equation (2) to 

obtain new values for the total strains e^,, £g , and e z . 
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In some cases it is possible to expedite the calculation of plastic 
strains for known values of total strains by preparing charts in advance 
of the calculation. Use of such charts will be illustrated in the exam- 
ples for the thin circular shell. 

The question arises as to whether the process is always convergent 
to the correct solution , or whether it is possible at seme point for suc- 
cessive solutions to become worse than earlier ones, thus they could lead 
intimately to meaningless results. In all the cases treated in this re- 
port it can be shown that the functions appearing in the integral equa- 
tions satisfy the conditions necessary for the convergence of the process 
(ref. 2); hence, the correct solution must result if a sufficiently large 
number of iterations are performed. However, for seme problems the rate 
of convergence may be very slew, and a large number of iterations may be 
required to obtain a solution of the desired accuracy. Two devices may 
be employed to expedite the convergence. It appears reasonable to assume 
that the closer the initially assumed solution is to the correct solution, 
the fewer the number of iterations that will be required before conver- 
gence will result. Hence, any knowledge or insight possessed by the in- 
vestigator should be used to estimate the zeroth approximation, rather 
t han obtaining it by assuming all plastic strains to be zero. Solution 
to related problems, solutions by other approximate or simpler methods, 
approximate measurements, physical intuition, or other expedients acces- 
sible to the investigator may all be used to good advantage. For example, 
same problems may be formulated in which the stresses will evidently 
c han ge little because of plastic flow, and the strains will have to assume 
whatever values are necessary to permit these stresses to be generated. 
(This is a contrasting case to that of thermally induced stress in which 
the strains govern and take on approximately their elastically c caput ed 
values, and the stresses adjust accordingly.) In such cases It may be 
better to start with an assumed initial stress distribution, compute the 
corresponding elastic strains from the stress -strain relation, and pro- 
ceed with successive integrations based on this initial strain 
distribution. 

Even when no insight into the correct solution is available, conver- 
gence may still be expedited by noting essentially the rate at which suc- 
cessive iterations change the strain distribution. A formal technique 
using this concept is illustrated in this report in connection with the 
rotating disk. It also should be noted that, if high-speed computing ma- 
chinery is used, the number of successive approximations required for 
convergence becomes of lesser importance. The method will now be illus- 
trated for four different problems . 
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EXAMPLES 
Thin Flat Plate 


As a first example , consider the simple uniaxial case of a thin in- 
finite plate of width 2c with a temperature distribution T(y) across 
the width. Under these conditions, the only nonzero stress is 
a = cr„(y) . As in the usual theory of Lending, it is assumed that plane 

sections remain plane. This requires that 

Sjj. = a + by ( 6 ) 

where a htiH b are constants to be determined. 

The stress -strain relation is 

Sc = 5 Sc + + Sep (7) 


V 


boundary conditions require that 

£ ‘ °J 

Combining equations (6), (7), and- (8) gives 

f E(a + by - oT - e^dy = 0 

and 

f E(a + by - aT - e^Jy dy = 0 
If E is constant, equations (9) give 


> 


6 dy 

xp 


—*£****£ 

j 


( 8 ) 


(9) 


( 10 ) 
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Also, in this case 


€ xp € ep 

a = cr 
e x 


Equations (6) can now be written as follows: 



As a specific example, let 

T 

E 


600 (y 2 - 1/3) + T q 
28X10 6 


c » 1 

a. = 9.5XL0 -6 


( 11 ) 


( 12 ) 


The stress-strain curve for the material 
of symmetry only half the plate between 
Equation (12) becomes 

e x - oT = e = -0.00570(y 2 - 


is given by figure 1. Because 
0 and 1 need be considered. 



(13) 


Equation (13) is now solved by successive approximations by using 
the stress-strain curve (fig. l) . For the zeroth approximation to the 
total strains it is assumed that the total strains do not change much 
from those computed "elastically" and, therefore, that e e p is every- 
where zero. The integral in equation (13) thus vanishes, and the strains 
are computed. This is the elastically computed strain distribution. 

With this strain distribution, a first approximation to the plastic 
strains e„ is obtained directly from the stress-strain curve (fig. l) . 

The integral in equation (13) is now evaluated, and a first approximation 
to the total strains is obtained. With this first approximation, a sec- 
ond approximation to the plastic strains is obtained from the stress- 
strain curve. The process is repeated until the desired convergence is 
obtained - The integral in equation (13 ) was evaluated simply by the 
trapezoi dal rule. More accurate evaluation, for example, by Simpson's 
rule, can be made if desired. 
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The computations for this problem are shown in table I, and the re- 
sults are plotted in figure 3. The stresses, which are not shown in 
table I, can be read directly from, the stress-strain curve once the strains 
are computed. These calculations show that the first approximation is 
sufficiently accurate and that the total strains are not much different 
from those computed elastically. This last result explains the fast con- 
vergence of this method for this problem. 

If the modulus of elastic ity E is not constant, no additional dif- 
ficulty is added. Solution of equation (9) gives, for a and b: 


pc pc 

i = A ^ J EoT dy + A 1 J Ee^ dy - 

p c pc 

A 2 J EccTy dy - Ag J Ee^ dy 

pc pc 

= -A 2 J EcdT dy - A^ J Ef^y dy + 


pc pc 

A 3 J_ c EoTydy + Aj^ Ee^y 


where A^_, Ag, and A3 are numbers which are determined once and for all 
for a particular problem from the known variation of E with temperature: 



i 
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By using these values for a and b the solution is obtained by succes- 
sive approximations as before. Of course, a different stress-strain 
curve must be used at every station and is dependent on the t emp erature 
of that station. 


Thin Circular Shell 


The second example considered is that of a thin circular shell with 
axial temperature gradient. The equations for the total stra ins , pre- 
sented in appendix C, are 


x R ' X 2 dx 2 


e z " 1 - V 


e e “ “ I 


( 6 x + e 0 ) + ^ - t “Hr + e 0p ) 



(Sep + ve ev )6z 




> (16) 


1 - 2v 


where w is the solution of 

d 4 v 


J 


d 2 ? 


. + 4w = -4EaT - Q 

dx 4 dx 2 


(17) 


and is given by 


w = c^ cos x cosh x + Cg cos x sinh x + c^ sin x cosh x + 


c^ sin x sinh x - 4R 


aT(!)G(x - |)d£ 


£ 


p(i) 


d 2 G(x - j) 


dx 


jr 

d| ’ Jo Q ^ )G ^ X " (18) 
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and 


G(x - g) = ^ 



g)cosh(x - - cos(x - g)sinh(x - g) 


1 



The term G(x - g) is the Green's function for this problem. Substitu- 
tion verifies that equation (18) is a solution of equation (17). 


The solution to the problem is now obtained by successive approxi- 
mations starting with the assumption that and e Qp are zero. A 

zeroth approximation to the total strains can thus be obtained from equa- 
tions (16) and (18) . With these values of total strain, first approxi- 
mations to the plastic strains and e Qp can be readily obtained 

as will be shown. These first approximations to the plastic strains are 
substituted in equations (16), (l8) , and (19), and first approximations 
to the total strains e^., e^, and e z are obtained. The process is re- 
peated as many times as necessary to give the desired convergence. For 
every successive approximation the constants c-^ to c^ appearing in 

equation (18) must be evaluated in order to satisfy the boundary 
conditions . 


To obtain the plastic strains and once the total strains 

e x , g . q , and e z have been determined for any iteration, the stress- 

strain curve and equations (3) and (4) are used. In order to facilitate 
the computation of e and Eg p , a cross plot is made of the stress- 

strain curve as is shown in figure 2 by means of equation (5) as pre- 
viously explained. The process for obtaining the plastic strains 
and 6g p from the total Btrains now becomes relatively simple. For a 
given set of strains e x , e^, and e , the equivalent total strain 
is computed by equation (3). For this value of e Q ^, e gp is obtained 
from figure 2, and and e 0 p are computed by equations (4). It 

should be noted that the curve in figure 2 is very close to a straight 
line. The equation of the "best fit" straight line was therefore used 
in some of the computations. Representation of the curve by a functional 
relation is particularly useful if a given problem is set up for auto- 
matic machine computation. 
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As an alternate procedure, the plastic strains e^ and can 

be obtained directly frcsn the total strains e^. and eg by means of a 
parametric family of curves as shown in figure 4. For a given pair of 
values of e x - a£E and - a,T, the plastic strains are read directly 

from these curves. Figure 4 was obtained from the stress-strain curve 
(fig. l) as explained in appendix D. This procedure avoids the neces- 
sity of ecuputing e z frcsn the last part of equations (16), e t from 

equation (3), and and e^ frcsn equations (4). However, obtaining 

a set of curves such as those in figure 5 involves a considerable amo unt 
of labor, and it is usually not worthwhile to make such a c har t unless 
several similar computations are to be made using the same stress-str ain 
curve. 


It is to be noted that the stresses can be computed at any step of 
the calculation by the stress-strain relations: 


0g = 




(eg - cdT 


■0p' 


(e x - aT - 


( 20 ) 


The solution will now be illustrated for a specific problem. Con- 
sider a thin circular shell with the following geometric and physical 
properties : 


L = 48 in. 
E = 12 in. 
H = 2 in. 


v = 0.3 


l 

E 


4 


e 2 h 2 


v 2 ) 


= 3.81 in. 


[3(1 

28X10® lb/sq in. 


a = 9.5XL0" 6 c y 1 
2 

T = 2.21x i (corresponds to 350° F rise from 
one end of shell to other) 


0 ^ x ^ 12.6 
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With the stress-strain curve of figure 1 and the following boundary 
conditions 


w(0) = w'(0) = v(l2.6) = w'(l 2.6) = 0 
the functions P and Q become 


( 21 ) 


P = 21 


•r 




( € xp + °* 3e 0p) z 42 


Q = 24 


s: 


> 


€ e P dz 


( 22 ) 


J 


Pram the first two boundary conditions. 


c^ = 0 


c 2 c 3 


Also, the first integral on the right side of equation (l8) becomes 
4K J* oT(g)G(x - |)d| => 0.000252x^ + 0.000252 sin x sinh x 
Equation (18) now becomes 

w(x) = Cg(cos x sinh x - sin x cosh x) + c 4 sin x sinh x - 


0.000252x 2 - I x (x) - I g (x) 


(23) 


where 


j i~X q(i)g(x " 




P(6) - d - G fe -- ^ d| 


(24) 


J 


For the zeroth approximation, it is assumed that e and and, 

therefore, and Ig are zero. The function w(x) is calculated from 
equation (23) with the constants eg and C4 determined from the last 
two boundary conditions (eq. (2l)). The strains eg, and e z 


are 
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then computed, frcm equations (16). First approximations to and 

are now obtained either directly from, figure 5 or by computing 
from, equation (3) , reading e e p from figure 2, and calculating 
and £gp frcm equations (4) . For this particular example figure 5 was 

used. For the two examples to be discussed subsequently, equations (3) 
and (4), in conjunction with figure 2 , were -used. With these values of 
and e^, P and Q are computed from equations (22), w from equa- 
tions (23) and (24), and new values of e JC , e Q , and e z from equations 
(16) . Hie process is repeated until convergence is obtained. 

Wherever derivatives of w are needed such as in the first part of 
equations (16) and in equation (21), the following relations are useful: 


aii r 

4 x 

* r= j 0 

Q(i) 

di„ 


P(§) 

d 2 i-. r 


^-1 

Q(l) 


dG(x - I) 
dx 






(25) 


d g G(x - g) 
dx 2 




dx 2 


= -4 



P(6)G(x - g)dg + P(x) 




and 


d 2 G(x - j) _ 


dx 


dG ^ X dx ^ = i sia ^ x ~ g)sinh(x - g) 

|sin(x - g)cosh(x - g) + cos(x - g)sinh(x 
d G(x - g) _ cos ( x _ |)cosh(x - g) 


*\ 


-] 


>(26) 
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The integrals in equations (24) and (25) were evaluated using the 
trapezoidal rule. Thus., in order to evaluate the integrals at a station 
x = 

^ fUjF&Cj- S)dt^(| f 0 B' J + f 1 P j . 1 + f 2 E ' J . 2 + ...+ fj» 0 ) 

(27) 

The results for this problem are shown in figure 5. An abbreviated 
calculation setup for one of the successive approximations is shown in 
table l(b). As many as seven successive approximations were carried out 
as shown in the figures , and the differences between the seventh and the 
fourth approximations are very small. From an engineering viewpoint the 
first approximation is actually sufficient. Again, the total strains 
did not change very much, which explains the relatively rapid convergence. 


Long Solid Cylinder 

The plane strain problem of a long solid cylinder with a radial tem- 
perature distribution is considered next. Appendix E shows that if the 
modulus of elasticity E is assumed constant, the total strains are 
given by 
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Equations (28) are now solved by successive approximations as for 

- the previous examples. The zeroth approximation to the total strains is 

obtained from equations (28) by assuming that e and e a are zero. 

jrp yp 

The total equivalent strain is then computed by equation (3), the equiv- 
alent plastic strain is read from. figure 2 , and first approximations to 
and £g p are obtained from. equations (4). These values of 
&£id Ggp are substituted into equations (28), and new approximations 
are obtained for e 0 , and e z . The process is then repeated as many 

gj times as necessary to obtain the desired degree of convergence. The 

lo stresses can be computed at any time from the general stress-strain 

^ relations: 


A ( e r + e 9 + £ z 
A(s r + e Q + e z 

A(e r + e 0 + e z 


- 3cc,T) + 2G(s r 

- 3ctT) + 2G(e 0 

- 3oT) + 2G(e z 


- ctT 




" 031 ~ € 0p) 

- aT + + e 0p ) 


(29) 


The above calculations have been carried out for a 1- inch- radius 
cylinder with a temperature gradient as shown in figure 6 and by using 
the stress-strain curve of figure 1. The ccmputat ions are shown in 
table 1(c) for one iteration, and the results are plotted in figure 7. 
Little difference occurred between the fourth and fifth approximations. 


Rotating Disk with Temperature Gradient 

As a final example the plane stress problem of a parallel- sided ro- 
tating disk with a radial temperature gradient and a constant value of 
E will be considered. In this problem the strains are considerably 
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larger than those in the previous examples . The total strains for this 
case, as shown in appendix E, are given by 


1 - pfl^r^ 

6 0 ~ " E 8 + 


i + v r 1 

Jo 


^1 




~rp ~ € 0p 


dr + 


(e„^ + e, )r dr + 
' rp 0p' 


€ r = " € 0 " “ '"I - " ^ + ^ + v)aT + e rp + V6 0p + 


(i-v) r hsL^Zez 
do 


dr + C, 
r 3 


e i — 

z 1 - v 


(«r + e 0^ + r^T: ^ - V T C e ^n + 6 flJ 


C, = 2(1 


- V) [$ ♦ 5£JL p** ♦ 


v 1 - V v rp T 0p^ 

R nR 


r 

v-^0 


*R 


2R' 


(€r * + e <* )r 


■] 


aTr dr - f / dr + 




(30) 

The solution to this problem is obtained by successive approximation 
exactly as in the previous example for the solid cylinder. By starting 
with assumed values of and Sgp equal to zero, e r , eg, and e z 

are confuted from equations (30), e et from equation (3), e g p from fig- 
ure 2, and e^ and e^ from equations (4). New values of e r , eg, 
and e_ are now obtained from equations (30), and the process is 
repeated. 

A solution was obtained in this manner for a 10- inch- diameter disk 

2 

with a temperature gradient as shown in figure 8. The value of pen was 
t ak en as 1500, and the stress-strain curve of figure 1 was used. The 
computations for one iteration are shown in table l(d) , and the results 
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are plotted in figure 9 , In this problem, the strains are relatively 
large, with the maximum, equivalent strain close to 1 percent. A straight- 
forward application of this method therefore requires approximately 40 
iterations in order to obtain accurate results. However, the convergence 
can be greatly increased by performing three or four iterations, taking 
the differences between successive iterations for the various strains, 
and extrapolating to a zero difference as shown in figure 10. A root- 
mean-square line is drawn, and the intercept at zero Ae is obtained. 

This furnishes anew starting estimate. Three or four more successive 
approximations are carried out, and another simi l ar extrapolation is made. 
This technique reduced the number of successive approximations for this 
problem from, about 40 to about 12. 

This same problem was solved by the trial-and-error method of ref- 
erence 3. The results obtained were almost identical to those obtained 
herein as can be seen in figure 9. 


DISCUSSION 

The speed of convergence of this method depends primarily on two 
factors: the amount of plastic flow occurring, and the number of stations 

taken in the plastic region. For small plastic strains convergence will 
be relatively fast. Similarly, for a small number of stations in the 
plastic region relatively few iterations are needed. If the number of 
stations is increased, more iterations are needed for convergence to 
occur. Thus, for the case of the thin shell of the second example, dou- 
bling the number of stations approximately doubles the number of itera- 
tions required for convergence. This is due to the fact that a change in 
the approximation at one station changes the values at all the other sta- 
tions, anti the more stations there are the longer it takes for all the 
stations to converge. 

Increasing the number of stations therefore increases the labor re- 
quired for two reasons: the time per iteration goes up, and the number 

of iterations required increases. Of course, the greater the number of 
stations used, the greater is the final accuracy attainable. However, 
good accuracy can be obtained without increasing the labor greatly by 
following one or both of the following techniques . A calculation is made 
with a Rmai 1 number of stations to obtain an approximate plastic strain 
dis tribution - This strain distribution is then used as a first approxi- 
mation in a new calculation using more stations. Also, an extrapolation 
technique such as shown in the rotating disk example can be used to speed 
up convergence greatly. It should be noted that doubling the number of 
stations for the thin- shell problem increased the labor greatly, but the 
maximum stress was changed by only about 2000 pou nds per square inch. 
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An interesting possible application of the results of the method 
presented herein may be worth further investigation. It has been found 
that, for thermal stress problems without additional loads, the total 
strains do not change very much because of the plastic deformation com- 
pared to the elastically computed strains and that the first approxi- 
mation to the stresses is thus usually fairly good. This leads to the 
possibility of determining s emiemp ir ically the plastic thermal stresses 
in complicated structures for which even the elastic stress distribution 
cannot be calculated and for which strain measur eme nts cannot readily 
be made in the plastic region and at high temperatures . A model of the 
structure can be constructed, and temperatures and temperature gradients 
can be simulated on a proportionally reduced scale so that no part of 
the model flows plastically. The total elastic strains under these con- 
ditions can usually be readily measured by means of strain gages. The 
elastic strains can then be extrapolated by simple proportion to those 
that would exist at the higher temperatures and gradients actually exist- 
ing in the structure. By assuming that the total strains are then equal 
to the strains that would exist if the material remained elastic, the 
plastic strains are computed from equations (3) and (4) and the stress- 
strain curve, and an estimate of the stresses is obtained from the stress 
strain relations such as equations (29). From the examples presented 
herein, it would seem that the stresses computed in this manner should 
be accurate enough for many engineering applications . 

In many practical problems it is necessary to take into account pre- 
vious plant ic flow that may have taken place. Thus, for example, in a 
thermal shock experiment, plastic flow may start at some time during the 
quenching process, and the material may continue to flow plastically as 
the process continues. A solution must therefore be obtained at various 
time increments from the start of the quench, and for each time interval 
the plastic flow that has already occurred up to that time must be taken 
into account , Treatment of this case is described in detail in appendix 
G. 


Although the method has been presented for four specific thermal 
stress problems, it is apparent that it is general in nature and can be 
applied to a large variety of problems for which the solutions of the 
elasticity equations are available. No implication is intended, however, 
that this method is necessarily more accurate or faster than other meth- 
ods that might be used for specific problems. Thus, the example of the 
rotating disk with temperature gradient can be treated more rapidly by 
the method of reference 3. The method of this report does, however, pro- 
vide a uniform simple approach that can be used for many different types 
of problems. It is not necessary, therefore, to develop special methods 
and techniques to handle different types of problems. 

Finally, it should be pointed out that, although for uniformity and 
simplicity the method has been set up by using the equations for total 
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strains (e.g., eg. (2)), it may be desirable in some cases to deal with 
the equations for stress which can be put in a similar form.. Upon deter- 
mining cr r , Cq, and a z for e^ = = 0 (i.e., the elastic stresses 

the elastic total strains can be computed by Hooke’s Law. For these 
total strains the plastic components are evaluated as described earlier, 
and a first iteration for stress is obtained. Subsequent iterations may 
follow the same procedure by using as the plastic strains the values de- 
termined from the previous iteration. In those cases where it is sus- 
pected that the elastic stress distribution is likely to be less affected 
by the plastic flow than the strain distribution, the plastic strains are 
determined directly from the stresses and the stress-strain relations as 
given in appendix B. 


SUMMARY CF RESULTS 

A method has been presented for solving plastic deformation problems 
in the elastic -plastic range. The method, one of successive approxima- 
tions, is illustrated by four examples which included a flat plate, a 
thin shell, a solid cylinder, and a rotating disk. It was found that for 
thermal stress problems accurate answers could be obtained with relatively 
few successive approximations. A technique for speeding up convergence is 
also shown. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 20, 1957 
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APPENDIX A 


Ai,Ag,A 3 

a,b 

C 1 ,C '2 ,C 3 ,C 4 


c l' c 2 ,c 3 jC 4 

E 

G 

G(x) 

E 

h 

I 1' I 2 

K i; K 2 

L ' 

l 

P(x) 

q(x) 

R 

r 

T 


SYMBOLS 

constants 

constants 

Integration constants 
half width of thin plate 
integration constants 
modulus of elasticity 

E 

Lame's constant , '+ ~ v) ' 

function 

thickness of thin shell 
thickness of rotating disk 
integrals 

constants (eqs. (B6)) 
length of thin shell 

characteristic length of thin shell, A I — — 

\[sU - v 2 ) 

function 


T, 


0 


u 


function 

mean radius of shell, or radius of solid cylinder 

radial distance to arbitrary point in solid cylinder 
or rotating disk 

temperature above arbitrary zero 
arbitrary constant temperature 

axial displacement of point on middle surface of shell 
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w 


x 


y 

z 


a 

6 


radial displacement of point cel middie surface of shell, 
positive inward 

ratio of axial coordinate of shell to characteristic 
length, or axial coor din ate of thin plate 

distance along width of thin plate 

radial coordinate of thin shell measured from middle 
surface, positive inward, or axial coordinate of 
long solid cylinder 

linear coefficient of thermal expansion 

conventional strain in tensile test 


e 

ep 

e et 

S:' e 9’ e z ,e T 
e xp ’ € 0p ’ e zp’ e rp 
9 


equivalent plastic strain 
equivalent total strain 

strains in x-, 9 - , z-, and r-directions, respectively 
plastic parts of e x , e 0 , e z , and e r , respectively 

tangential coordinate 


A 

v 

i 

P 

ff e 

CT x jO 0' °r> °z 
a> 


Lam/'s constant, 2 v) 

Poisson's ratio 
integration variable 
density 

equivalent stress 

normal stresses in x-, 0-, r-, and z-directions , 
respectively 

rotational speed of disk 
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APPENDIX B 


CALCULAnUCN OF PLASTIC STRAINS 

The deformation theory of plasticity is used with the three usual 
Assumptions that the directions of the principal strains coincide with 
the directions of the principal stresses, that the ratios of the prin- 
cipal shear strains are equal to the ratios of the principal shear 
stresses, and that the volume remains constant in the plastic range. 
These assumptions imply 



into the first part of equations (Bl) , it can also be shown that 



where 

% = K 2 + ^-±-2 ( B 4) 
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Define 



Then, by squaring and adding the equations in (Bl) and (B3) it readily 
follows that 



Hence, by the relation between Kj and Kg in equation (B4) 


€ et = 


2(1 + v) °e 


The plastic strains can be determined in terms of the total strains by 
dividing equation (B3) by equations (Bl) and applying equations (B6): 



Solving equation (B7) and the incompressibility relation in equations 
(Bl) results in equations (4). 
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APPENDIX C 


EQUATIONS FOR THIN CIRCULAR SHELL 

The equilibrium equations for a thin circular cylindrical 
given in reference 4: 

N x = 0 


R Ac 


Z 2 dx 2 6 


+ N 0 = 0 


■where 


**S* 


N„ 


H/2 

H/2 

H/2 




°x 42 


\fn/ 2 ° e 32 


pH/2 

a > 

The stress-strain relations, including the plastic strains 

"x = ^7-^ [s - aT - Sep + v C e 0 - ^ - e 0p^ 

a 9 " ” 7 |^0 -an - £gp + v(e x -an - exp)j 


,z dz 


J 


The strain displacement relations are 

du 


=2. t 
d v 

X “ l dx " z 2 ^2 


€„ = 


> 


W 

e 6 = ” R 


J 


shell are 
(Cl) 


(C2) 


, are 
( 03 ) 


( 04 ) 
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ii 

■ 

O 

n 


Substituting equations (C4) into (C3) gives 




du d 2 y 

i dx " ,2 , , 
l dx 




2 z ' 01 ' Sp + v (' j ‘ ‘ cdr ) 


aT - e 0p + v 


f du d 2 v 
V Z to " Z 2 dx 2 


z - a,T - 


Sq) 


>CC5) 


J 


Frcm the first equation of equations (Cl) and (C2), 

dU V | + (1 + v)al + j f/ <Sp + ''V dz 


1 dx 


(C6) 


H 9 = -EH 


(i * i r 

V R E J-H/2 


e 0p 42 + aT 


) 






EH 3 


d 2 w 

2 ,2 .2 


E 


vH/2 


> (C7) 


12(1 - vr l dx 1 - v 


C £ xp + ve 6p) z 32 


y-H/2 j 

Substituting equation (C6) into the first of equations (C4) gives 

\H/2 
-H/2 

Also, from the second of equations (Cl), 


= v g + d + v)aT + g ( 

U-H/2 


£ *p * * - pfel 


z (C8) 


i/2 


d 4 w ^ 12 Z 2 d 2 / , v J 4H 

dx 4 ' H 3 dx 2 ^ Ve0 l ,)z " H 


To obtain e z , the third stress -strain relation is used: 



:/ 2 


'-H/2 


= 1 j^z - v(a x + ff 0 ) + € zp + off] 


€ ep 32 


(C9) 


(CIO) 


Substituting a z = 0 and equations (C3) into equation (CIO) gives the 
third part of equations (16). 
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APPENDIX D 


PLASTIC STRAIN CHARTS 

In order to obtain the chart shown in figure 4, the stress-strain 
relations with a_ = 0 are written as follows: 

Zi 


X 


- = B (a x - vo e> + e xp 


- or = I (c ( 


E ^ 9 


vff x) + 6 0p 


(Dl) 


Also, 


3 € ep ) + € ep) 




^ (2e 9P + V, 


> 


0 


Substituting equations (D2) into (Dl) gives 


aT = 


e - aT = 
0 


1 + 


2 (2 - y) g e 


-ep 


1 + 3 


2 (2 - v) ff e 1 
E 


ep 


Sep 


2 (l - 2v) °e 




E 


e e 0p 
ep 


e 0p + 3 


2 (l - 2v) a e 

E €_ Sep 


-ep 




(D2) 


> (D3) 


With the above equations and the stress-strain curve for the mate- 
rial, a two-parameter family of curves can he plotted giving the total 
strains for any pair of plastic strains and Thus, 

(1) An arbitrary convenient value is chosen for e^. 

(2) A series of values are chosen for . For each of these values 

(a) Compute e or , from equations (B5) . 

(b) Read from the stress-strain curve. 

(c) Compute ^ - cdT and e 0 - aT from equations (D3). Thus, 

one curve of the family is obtained. 

(3) To obtain the other curves, new values are chosen for e^p, 
the process is repeated in each case. 


and 
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The limiting curve of zero plastic strain is an ellipse about the 
origin as shown in figure 4. Any point inside this ellipse corresponds 
to zero plastic strain. 
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APPENDIX E 


LONG SOLID CYLINDER 


Consider an infinitely long circular cylinder which has a radial 
temperature distribution T(r) . The equilibrium equation is 


da a - a - - 

U + « . o 

dr r 


(El) 


The compatibility relation is 


id€c 

e a = r — 
9 dr 


(E2) 


Substituting equations (29) into (El), assuming E constant, and elim- 
inating e r by use of equation (E2) give 


_d_ 

dr 


1 d_ 

r dr 




(*■%) 


1 + v dT.l-2vld 


a 


dr 1 - v 


) 1 - 2V 9p 

r dr ^ rp J 1 - v r ( ‘ E3 ' ) 


Integrating equation (E3) results in 

vr 


"0 


1 ' V r 2 Jo 1 - v r 2 Jo 


e^r dr + 
rp 


1 - 


fif’I 


T 6 5P -fg£ dr dr + C 1 + ^| (E4) 

r r 2 


For a solid cylinder Cg must equal zero. Also, equation (E2) can be 
rewritten as 


1 d r 2 \ 

e_ = + - — (r e Q ) 

r 0 r dr 9 


(E5) 


By using equation (E3), equation (E5) can be written as: 

*r 


6 r “ - e 9 + i _ v 




/ 


fg£JLfg£ 


dr + 2C-, 


(E6) 
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To determine use is made of the fact that e_ is a constant and 

that z ■ 

R 

c z r dr = 0 (E7 ) 

Substituting the last of equations (29) into (E7) and using equations 
(E4) and (E6) enables e z to be determined. To obtain C the surface 

boundary condition cr r (R) =0 is used. Substituting the first of equa- 
tions (29) into this relation enables to be calculated. The final 

results are given in equations (28) . 
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APPENDIX P 


ROTATING DISK 


The derivation of the equations for the rotating disk with a radial 
temperature distribution is very similar to the derivations of the equa- 
tions for the long solid cylinder (appendix E) . Eor the disk problem, 
the axial stress 0 Z is assumed to be zero, and the axial strain e z 
is no longer a constant. The equilibrium equation for this case is 


(hro r ) - ho fl + paAr 2 = 0 


(FI) 


The compatibility equation is the same as equation (E2), and the stresB- 
strain relations with cr z equal to zero become 


*r = 


- e r 

= — T V - V T 

x . v 2 [r 0 ^ rp 0p 


- 


E 


1 - v" 


+ ve fl - (e^ + ve^J - (l + v)arj 
|^0 + V€ r - ( € 0p + V€ rp) - ^ + v > a ^] 


(F2) 


For a parallel- sided disk with E assumed constant, the solution 
is readily obtained. Substituting equations (F2) into (Fl) and using 
equation (E2) result in 


_d_ 1 _d_ 

dr r dr 


(r 2 e 0 ) 


. . + (1 + ,) iJgl ♦ 

d , ^ , e rp " e 0p 

d? (e rp + ve ep } + ^ - v ) r 


(F3) 


and, after integrating, 


(l - v 2 ) pcn 2 r 2 1 + v 
= - E + ~ ^2" 


aTr dr + 


1 - v 


f 

'Jn 


~rp ~ 0p 


2r v/o 


v C 3 C 4 

e rp + W 41 + T + ^ 


dr + 


(F4) 
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where for a solid disk most vanish. To obtain equation (F4) 

is substituted into equation (E2) resulting in 


£ r " - e 0 


- T ^ - (1 **)"* « Sp ♦ »Vp + 


(1 - v) 


X 


ZblzJSs . + c 

r 3 


(F5) 


The term. e z can now be computed from equations (B2) with c z taken as 
zero. Thus, 

v r \ 1-2 v r 

€ z - - 1 _ v ^ e r + 6 0 ^ “ i _ v 


e iT + 


w + Hr 01 


Ore) 


The constant C 3 is evaluated from the known rim loading. If the rim 

stress due to the rim loading is a r (R), the first of equations (F2) 
becomes 


<J r ( R ) = "" 2 S + V€ 0 + (Sp + ve 0p) ‘ C 1 + v ) aT 

1 “ v L -Jr=H 


(F7) 


Substituting equations (F4) and (F5) into equation (F7) enables C 3 to 
be calculated. The final equations are given in equations (30) . 
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APPENDIX G 


EFFECT OF PREVIOUS PLASTIC STRAIN 

If the body under consideration has undergone previous plastic 
strain such as that during a thermal shock in which the stresses and the 
plastic deformations vary with time, a complete series of solutions must 
be carried out for successive time intervals. The equations for the 
strains remain the same except that each plastic strain term such as 

is replaced by + Ae^, where le^ is the total accumulated 

plastic strain up to the previous time and is the additional plas- 

tic strain during the time interval under consideration. These strains 
Ze^ and 2 e gp are known from the previous calculations, and the total 

strain equations are solved by successive approximation to obtain the 
change in plastic strains and Ae„ as well as the total strains 

rp op 

after the new time interval . 

In carrying out this type of calculation for successive time inter- 
vals, a time may eventually be reached when some point in the body begins 
to unload, that is, starts decreasing. When this time is reached, 
no additional plastic flow will take place at this point, and unloading 
will proceed along an elastic line such as CC 1 in figure 1. At this 
station then, the plastic strains are assumed to be zero from this time 
on. The successive approximations are continued in the usual manner until 
all points in the cylinder begin to unload or equilibrium conditions are 
reached . 

This type of procedure of adding the plastic flow occurring after 
each time increment to the previously accumulated plastic flow is equiv- 
alent to the assumption that, as the load and temperature change, the 
stress position on the new stress-strain curve would be the same as if 
a test specimen were loaded above the yield point, the load removed, the 
temperature changed, and a new load applied. This assumption is illus- 
trated in figure 11 in which point A represents a loading at the first 
temperature conditions; the dotted line AB represents the unloading path; 
the curve BCD shows the stress -strain curve at the new temperature; and 
point C gives the new stress position. The total strain at this point C 
is given by the sum of three strains: the residual strain caused by the 

first loading, the elastic part of the strain caused by the second load- 
ing, and the plastic strain caused by the second loading. 

When the foregoing procedure is applied, the curve BCD must, of 
course, represent the true stress-strain curve at the new temperature of 
a material that has already been subjected to the plastic cycle GAB. In 
general, this new stress-strain curve is different from the stress-strain 
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curve at the given temperature of a material that has not been subjected 
to plastic flow. However, unless data are available, it may be necessary 
to assume that the curve BCD is the stress-strain curve at the given 
temperature of a specimen of virgin material. Results obtained in this 
way, however, should he treated cautiously since this effect may he 
appreciable. 
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TABLE X. - PLASTIC FLCW CALCULATIONS 
(a) Flat plate 


Station 

7 

®ep 

( Bq . ‘( 13 )) 

'IJP 

(Fig. 1) 

( Bq . ( 13 )) 

e ep 

(Fig. 1) 

(Eq. (13)) 

/ 6e]p 
(Fig. 1) 

( Eq . ( 13 )) 

%p 

(Fig. l) 

( Eq . S ( l3 )) 

e ep 

(Fig. 1) 

1 

0 

0 

1 . 9XL0 -5 

0 . 7XKT 3 

1 . 7X10- 3 

0 . 5XL0" 3 

1 . 6X10-3 

0 . 4X10- 3 

1 . 6X10- 3 

0 . 4X10- 3 

1 . 6X10-2 

0 . 4x10-2 

2 

.1 

0 

l.a 

.6 

1.6 

.4 

1.5 

.4 

1.5 

.4 

1.5 

.4 

3 

.2 

0 

1.7 

.5 

1.5 

.4 

1.4 

.3 

1,4 

.3 

1.4 

.3 

4 

.3 

0 

1.4 

.3 

1.2 

.2 

1.1 

.1 

1.1 

.1 

1.1 

.1 

5 

.4 

0 

1.0 

0 

.6 

0 

.7 

1 0 

,7 

0 

.7 

0 

6 

.5 

0 

.5 

0 

.3 

0 

.2 

0 

.2 

0 

.2 

0 

7 

.6 

0 

-.2 

0 

-.4 

0 

-.4 

0 

-.5 

0 

-.5 

0 

a 

.7 

0 

-.9 

0 

- 1.1 

-.1 

- 1.2 

-.2 

- 1,2 

-.2 

- 1.2 

-.2 

9 

.6 

0 

- 1.7 

-.5 

- 1.9 

-.7 

- 2.0 

-.8 

- 2.0 

-.8 

- 2.0 

-.8 

10 

.9 

0 

- 2.7 

- 1.5 

- 2.9 

- 1.6 

- 3.0 

- 1.7 

- 3.0 

- 1.7 

- 3.0 

- 1.7 

11 

1.0 

0 

- 3.8 

- 2.5 

- 4.0 

- 2.7 

- 4.0 

- 2.7 

-4.1 

- 2.8 

- 4.1 

- 2.8 


(b) Thin shell ■ 7CS-150T6 . 
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TABLE I. - Concluded. PLASTIC PLOW CALCULATIONS 


(c) Long solid oyllnder 


Station 

r 

aT 

e rp 

(Previ- 

ous 

approxi- 

mation) 

60 P 

(Previ- 
ous " 
approx! 
mation) 

e r 

{Eq. (28)) 

B e 

(Eq. (28)) 

B z 

(Eq. (28)) 

(Eq. S (3)) 

. e P 

(Fig. 

2) 

6 rp 

(Eq. 

(4)) 

e 6p 

(Eq. 

(4)) 

1 

0 

9.50XLCT 3 

0 

0 

9.457X10 -3 

9.457xlO“ 3 

8.912X10“ 3 

0.363XX0" 3 

o>acr 3 

0XL0 -3 

OXLCT 3 

2 

.75 

9.50 

0 

0 

9.462 

9.452 

0 

.363 

0 

0 

0 

3 

.80 

9.29 

0 

0 

9.085 

9.439 

0 

.310 

0 

0 

0 

4 

.85 

8.81 

0 

0 

8.240 

9.393 

0 

.668 

0 

0 

0 

5 

.90 

7.68 

0 

0 

6.236 

9.299 

0 

1.926 

0.91 

-.907 . 

0.545 

6 

,95 

7.06 

0 

0 

5,274 

9_1 Od 

0 

2,494 

1,40 

-1.397 

. 754 

7 

1.000 

6.83 

0 

0 

5.045 

8.911 

0 

2.578 

1.48 

4.477 

.738 


(d) Rotating dlak 


705-15078 


Dbti.LJ.UII 


— tn 
Ul 

! rp 

( Previous 
approxi- 
mation ) 

. ®6P 

(Previous 
approxi- 
mation ) 

(Eq. (30)) 

*9 

(Hq. (30)) 

(Eq. (30)) 

(Bq. B (3)) 

s ep 

(Pig. 2) 

®rp 

(Eq- (*)) 

r. _ 

-ep 

(Eq. (4)) 

1 

0 

0.95X1CT 3 

1.567x10"® 

1 .567X10"® 

3. 757x1 CT® 

3.757x10"® 

-3.247X10"® 

4.669X10"® 

3.441x10”® 

1.720X10"® 

1.720x10"® 

£ 

.5 

1.045 

1.515 

1.517 

3.780 

3.764 

-3.033 

4.543 

5.317 

1.657 

1.660 

3 

1.0 

1.53 

1.459 

1.361 

3.983 

3.845 

-2.490 

4.270 

3.048 

1.573 

1.474 

4 

1.6 

1.71 

1.376 

1.144 * 


3.942 

-1.776 

3.919 

2.721 

1.466 

1.252 

5 

2.0 

2.47 

1.256 

.7269 

4.862 

4.104 

-.3885 

3.277 

2.124 

1.298 

.8067 

6 

2.5 

3.80 

1.089 

.1111 


4.381 

1.913 

2.389 

1.304 

1.047 

.1508 

7 

3.0 

5.416 

1.070 

-.428 a 

7.614 

4.797 

4.371 

2.035 

.9972 

.9899 

-.3908 

8 

3.5 

7.315 

1.410 

-1.175 

9.773 

5.362 

6.964 

2.578 

1.469 

1.371 

-1.142 

9 

4.0 

9.5 

1.904 

-2.298 

12.28 

6.076 

10.00 

3.623 

2.44S 

1.908 

-2.279 

10 

4.6 

11.88 

2.290 

-3.637 

14.83 

6.911 

13.51 

4.696 

3.664 

2.303 

-3.620 

— 11 

5.0 

14.25 

2.498 

-4.996 

17.18 

7. 824 

17.18 

6.235 

4.982 

2.491 

-4.902 


W 
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STRAIN, € / CS- 15125 


Figure 1. - Typical stress-strain curve for 18-8 stainless 
steel. 
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( a ) Strain /CS-15086/ 


Figure 3. - Variation of stress and strain along 
flat plate for successive approximations. Tem- 
perature given by T = 600(y2 _ i) + Tq. 
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(b) Stress. /CS- 15088/ 


Figure 3. - Concluded. Variation of stress and strain along 
flat plate for successive approximations. Temperature given 

by T = SOOCy 2 - i) + T 0 . 
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Figure 4 . - Plastic strain chart based on Btress-strain curve of figure 1 . 
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(b) Tangential stress on outside surface of shell. 


Figure 5. - Continued. Variation of stresses and strains on shell with 
axial distance measured in terras of characteristic length. Tempera- 
ture distribution, T «« 2.21x^. 
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(c) Axial strain on outside surface of thin shell. 


Figure 5. - Continued. Variation of stresses and strains on shell with 
axial distance measured in terms of characteristic length. Temperature 
distribution, T = 2.21x 2 . 
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(d) Tangential strain on outside surface of shell. 

Figure 5. - Concluded. Variation of stresses and strains on shell vith axial dis- 
tance measured in terms of characteristic length. Temperature distribution, 

T = 2.21x 2 . 
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Figure 6. - Radial temperature distribution in long solid 
cylinder . 
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(a) Strain. 

Figure 7. - Strain and. stress distributions in long solid 
cylinder with radial temperature gradient. 
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(a) Strain. 


Figure 9. - Strain and stress distributions in rotating disk with 
temperature gradient, par = 1500. 
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(b) Stress. 


/CS-13092/ 


Figure 9. - Concluded. Strain and stress distributions in rotating 
disk with temperature gradient, pee 2 = 1500. 
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Figure 11. - Uniaxial Btress-strain curves showing 
components of strain when plastic flow occurs a 
second time. 
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